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We present a theory for the description of energy relaxation in a condensate of bosonic particles, 
taking the microcavity exciton-polariton system as an example. The approach is based on a Monte 
Carlo type treatment of the phonon field and a microscopic Hamiltonian for exciton-phonon inter- 
actions. This approach, which also fully accounts for the Coulomb (exciton-exciton) interactions, 
pumping and particle lifetime, offers a complete description of experiments in which Bose-Einstein 
condensates are seen to relax their energy as they propagate in real space and time. 
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Introduction. — Bose-Einstein condensates form when 
multiple bosonic particles relax their energy, through in- 
teraction with other particles (e.g., phonons), to collect 
into a low energy state. Examples in solid-state sys- 
tems include: condensates of magnons [l], in quantum 
spin gases; indirect excitons 0] in coupled quantum wells; 
and exciton-polaritons in semiconductor microcavi- 
ties. While Bose-Einstein condensation is conventionally 
thought of as a macroscopic occupation of the ground 
state in thermal equilibrium, an important issue is that 
thermal equilibrium in solid state systems is never per- 
fectly achieved since the particles have finite lifetimes. 
The departure from thermal equilibrium is particularly 
pronounced in exciton-polariton systems, where a short 
lifetime of particles may make the system highly non- 
equilibrium 0- Nevertheless, spontaneous spatial coher- 
ence has been observed Q , althourfi not necessarily in 
the ground state of the system [1, Q. The theoretical 
description of such condensates thus requires a kinetic 
approach. 

In excitonic systems (exciton or exciton-polariton con- 
densates) scattering with acoustic phonons offers a 
mechanism of energy relaxation. It may not offer com- 
plete thermalization, but the effect of energy relaxation 
is clearly seen in experiments where a potential gradient 
is present. The latter can be created by application of 
stress m, optically [lH or structurally engineered (l^ . 
In this situation, the processes of energy relaxation and 
condensate propagation down the gradient are closely 
connected. Aside offering a clear demonstration of en- 
ergy relaxation, propagating condensates have recently 
demonstrated potential as optoelectronic transistors in 
exciton Q and exciton-polariton systems [isj . 

Theoretically, the dynamics of spatially homogeneous 
systems have been described accounting for the exciton- 
phonon interaction by means of a system of semiclassical 
Boltzmann equations This approach, however, 

has serious drawbacks. First, the corresponding formal- 
ism is based on the assumption of full incoherence in the 
system, thus quantum states are supposed to be com- 
pletely uncorrelated. However, as long as the condensate 
has been formed, this is no longer valid. As a conse- 



quence, a number of coherent phenomena such as the 
onset of superfluidity bistability [20| and hysteresis 
cannot be described. Second, the Boltzmann equations 
can provide us with information about the occupation 
numbers in reciprocal (A:-space) only, whereas the real 
space (x-space) behavior remains obscure. 

On the other hand, if the processes of decoherence 
arc fully neglected, the state of interacting particles can 
be treated as a classical field described by the Gross- 
Pitaevskii equation Ll9|, which can be modified for in- 
coherent pumping [7|, |2l|. Such an approach has been 
successful for the description of a variety of recent exper- 
iments [2^ . [2^ in semiconductor microcavities, includ- 
ing, for example, experiments on the dynamics of vor- 
tices I22II. spatial pattern formation 23|, |2j| and spin tex- 
tures [25|, l26|. However, the Gross-Pitaevskii equation 
conserves the particle energy and thus does not account 
for phonon-assisted scattering. Recent models include 
energy relaxation in a phenomenological way within a 
classical stochastic field [27| or Gross-Pitaevskii type for- 
malism [H]. While these models give results in agree- 
ment with experimental data, they operate with un- 
known phenomenological parameters. 

In this Letter we introduce a microscopic theory for 
the description of energy relaxation in a coherent ex- 
citonic ensemble. The exciton/exciton-polariton field 
is coupled to a field representing phonons in the sys- 
tem, which is modeled using stochastic variables. For 
simplicity, we consider a resonant coherent injection of 
excitons/exciton-polaritons (models of non-resonant ex- 
citation, involving coupling to an exciton reservoir, have 
been considered elsewhere 0] and are compatible with 
our approach). We consider injection with zero in-plane 
momentum in a system containing a potential gradi- 
ent [13] ■ The potential gradients accelerate the parti- 
cles which undergo scattering with acoustic phonons as 
they propagate. The latter process leads to the energy 
dissipation and thermalization in the system. To give a 
complete description of the dynamics we fully account 
for exciton-exciton scattering and losses provided by fi- 
nite lifetime. 

We will use parameters corresponding to cav- 
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ity exciton-polaritons - the system for which non- 
equilibrium effects are most clearly pronounced - al- 
though it should be noted that our formalism is applica- 
ble to the modelling of other systems in which bosons re- 
lax their energy through interactions with an incoherent 
gas of other particles. Another example of application of 
our theory would be the system of indirect excito ns [3i . In 
the excitonic optoelectronic transistor (EXOT) (UH^l, 
electric fields introduce potential gradients and the con- 
trol of the electric fields modulates the fluxes of excitons. 
Similar to the exciton-polariton case, the relaxation down 
potential gradients should be mediated by phonons, the 
process that has not been addressed theoretically so far 
to the best of our knowledge. 

Theory. — For simplicity, we consider a ID system, as 
in the case of microwires where there has been a recent 
experimental focus on energy relaxation of propagating 
condensates ll|, |l3|. We note however that the for- 



malism is expected to be compatible with 2D systems 
as well. We introduce the quantum field operators for 
exciton-polaritons, connected with annihilation oper- 
ators in reciprocal (k-) space by the Fourier transforms 
{J- in short notation). 



where is the discretization length. 
The Hamiltonian of the system reads 



n 



q q,k 

The first two lines here correspond to coherent pro- 
cesses. Ek is the particle dispersion (which is non- 
parabolic for exciton-polaritons); the last term in the 
first line is the coherent pumping, where Vx is the in- 
tensity of pump, kp is the wavevector of the pump, cor- 
responding to the inclination of an incident laser beam, 
and hujp is the pumping energy; Vx is the potential pro- 
file in a;— space, a is a constant describing the strength 
of particle-particle interactions. 

The last line in the equation above corresponds to inco- 
herent processes. To model the interaction with acoustic 
phonons, we introduce the Frohlich Hamiltonian [lol . 31 1, 
where the phonons described by operators St, are con- 
sidered to be three dimensional. The phonon wavevector 
is q = CxQx + fiylx + ezQzj where e^, Cy and Cz are unit 
vectors: Cx is in the wire direction, is in the struc- 
ture growth direction. The phonon dispersion relation, 

is determined by the speed of 



^u^ql + q^ + qt, 
sound, u. is the exciton-phonon interaction strength. 



whose calculation can be found elsewhere |3ll - l35l | (see 
also the detailed derivation of the formalism in the "Sup- 
plementary Materials" appended to this Letter). 

Using the Heisenberg equations of motion for the op- 
erators, one can write the formal solution for the phonon 
field as following: 

kii) = MO)e-*'^f* 

'if 4+,. in it') e-^-^i'-'')dt'. (2) 

Remembering that phonons represent an incoherent 
thermal reservoir, we can replace the term 6g-(0)e~*'^9* by 
a stochastic classical variable hq^{i) (and a similar replace- 
ment can be made for the conjugate field). This repre- 
sents the analogue of the Markov approximation within 
the Langevin approach, when pho nons are assumed to 
have a randomly varying phase j32l | . The stochastic vari- 
ables are complex numbers with real and imaginary com- 
ponents normalized as follows. 



{h}{t)h,.{t'))^n^5,^^5{t-t')- 
{\~{t)b^,{t'))^{b}{t)h*^{t'))=Q, 



(3) 



where is the number of phonons in the state with 
wavevector q determined by the temperature of the sys- 
tem. 

The exciton-polariton field dynamics can then be de- 
termined solely by the exciton-polariton operators and 
the stochastic terms. Further, within the mean field 
approximation, the field operator can be treated 
as a classical variable for condensed exciton-polaritons, 
ipx = (^x) (with the Fourier image ipk)- Then, physi- 
cal observables are calculated over multiple realizations 
of the evolution dynamics with stochastic variables bkit)- 
The corresponding equation of motion reads (see the Sup- 
plementary Material for the details of the derivation): 



ih^ = [EkiJk + Sk{t)] 
at 

Vx+a \il)x\^ - 



^x + Vx^^^'^^r^'^^* 



(4) 



where we introduced phenomenologically the decay term 
—ihrijljxll to account for the radiative decay of parti- 
cles [l9[. The constant a describing polariton-polariton 
interactions can be estimated as dl: a w Ebal/{LyAx), 
where Ly is the lateral size of the microwire and Ax = 
Lx/N is the discretization unit. 

One can sec that interaction with phonons leads to the 
appcarence of two types of terms. First, one has a term 

Skit) = Y ^^+^^ (^) { f A. (^')^o(t - t')dt^ (5) 
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where Aq^{t) = J^k' "^I'+q^ (*) ^fe' (*)■ The term is pro- 
portional to the cube of the polariton field and does not 
directly include a stochastic term. It can be thus in- 
terpreted as the term corresponding to the emission of 
phonons by a condensate stimulated by polariton den- 
sity. The convolution integral is responsible for energy 
conservation. Note, that the function 

q,i.qz 

J \Gm^M^imdqydq. (6) 

is approximately independent of in the range of G 
(—10^,10*) m^^, and thus in our calculations we put 

iCgAt)-'Co{t). 

The stochastic functions Tn and 71* in the last line of 
Eq. ([M)) are defined by the correlators: 

QvQy 

(Or,i it')) = (r; (t)72 it')) = 0. (7) 

These thermal terms contain the phonon field and so are 
strongly temperature dependent. The proportionality of 
the thermal part to the first power of ipx in Eq. ([M)) 
means that the scattering processes proceed at a rate 
proportional to the first power in exciton-polariton den- 
sity. Consequently, this term corresponds to the absorb- 
tion of the phonons by the polariton ensemble and their 
emission which is stimulated by final state phonon oc- 
cupancy, but not exciton-polariton density. The latter 
processes are instead represented by the stimulated term 
(Eq. [5]) considered above. 

Results. — We considered an InGaAlAs alloy-based mi- 
crocavity and in computations used the following set of 
parameters: speed of sound u ~ 5370 m/s [31}, 7 = 1/18 
ps~^ [l3|- The function Vx ~ Vq ~ (3x is the exciton- 
polariton potential in space, composed of: the potential 
defining the walls of the exciton-polariton wire, Vq and a 
potential gradient /3 = 9 meV/mm. We consider a lin- 
ear potential gradient as in the experiment of Ref. [isj . 
Our formalism would however be compatible with any 
potential shape and could be applied to, for example, 
the parabolic [Hj or staircase [l^l potentials studied re- 
cently in experiments. The exciton-polariton dispersion 
was calculated using a two oscillator model with cavity 
photon effective mass 4 x 10~^ of the free electron mass, 
Rabi splitting 10 meV and exciton-photon detuning 2.5 
meV at zero in-plane wavevector. 

The main phenomenon we investigated was the re- 
laxation of energy of exciton-polaritons (thermalization) 
caused by the interaction with acoustic phonons. Exciton 
polaritons were introduced to the system by localized co- 
herent pump with energy coinciding with the bottom of 



lower polariton branch (zero detuning case) which guar- 
antees that we are outside bistable regime. 

Stochastic equation (|34p can be solved numerically. 
The results of the modeling are presented in Figs. [T][2^- 
f. Figure [1] illustrates the propagation of particles cre- 
ated by a short pulse along the ID wire. The pumping 
is switched on during the first 20 ps of the theoretical 
experiment, and then it is off and we observe the de- 
cay of intensity due to the finite particle lifetime. The 
quantity \tpx\'^ is depicted for different times: 20, 40, 60 
and 80 ps {t is the parameter). One can see that the 
wavepacket of exciton-polaritons propagates along the 
channel with time, disperses, and accumulates at the 
right hand side. Note, that reflection from the poten- 
tial jump at the end of the wire produces pronounced 
interference fringes clearly visible in recent experiments 

Figured] illustrates the energy relaxation in the system 
for three time intervals: 0-50 ps, 50-100 ps and 100-150 ps 
correspondingly. It is clearly seen that exciton-polaritons 
pumped at a; = 0, fc = (a,b) propagate along the poten- 
tial slope losing their energy (c,d) and accumulate on the 
right-hand side of the wire (e) thus condensing at fc = 
in reciprocal space (f). The plots (a,c,e) snapshot the 
x-space dynamics of the system, the fc-space dynamics is 
depicted in plots (b,d,f). The white dashed curves in the 
figures correspond to the potential profile (in x-space) 
and free polaritons dispersion (in fc-space). It should be 
noted, that the relaxation with phonons is an efficient 
process, thus particles rapidly move towards the bottom 
of the slope and accumulate there; and their concentra- 
tion near the pumping spot soon becomes lower than the 
concentration in the signal point. 

Discussion. — Comparing with previous theoretical 
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FIG. 1: Propagation of particles along the potential slope 
in the quantum wire due to phonon-assisted relaxation. The 
curves correspond to the particle concentration profile in x- 
space for different times: 20 (red), 40 (green), 60 (blue) and 
80 ps (magenta). The particles are introduced by a coherent 
Gaussian pump of the duration of 20 ps of the theoretical 
experiment. Particles propagate along the wire (green and 
blue curves), suffering losses caused by their finite lifetime 
and refiecting from the end of the wire. The interference 
of the incoming and reflected waves produces fringes clearly 
visible at 80 ps. 
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models aimed at description of energy relaxation in co- 
herent condensates, we would like to stress that all the 
parameters, including phonon scattering rates, are calcu- 
lated microscopically rather than being treated in phe- 
nomenological way. Furthermore our approach allows the 
possibility to treat large regions of direct and reciprocal 
space with reasonably low computational cost. This in 



in contrast to the theory developed in Ref. [28[ where 
the need to calculate an evolving particle spectrum in a 
coarse time window was extremely memory-demanding. 
In the approach of Ref. 127| the analysis was restricted 
to a single quantum state and the energy relaxation in- 
troduced phcnomenologically. Very recently, in Ref. [37| , 
an approach merging the Boltzmann equations with a 
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FIG. 2: Relaxation of energy of exciton-polaritons along the 
potential gradient in the regime of cw excitation in a quan- 
tum wire due to phonon-assisted processes in x- and fc-space 
for different time ranges: 0-50 ps (a,b), 50-100 ps (c,d) and 
100-150 ps (e,f). White curves show the model potential pro- 
file in the real space and free exciton-polariton dispersion in 
the reciprocal space used in calculations. Coherent pump- 
ing represents a Gaussian in x-space centered around x = Q. 
The inclination angle of the pump is zero, thus fcp = and 
fujjp = -E'(O). Polaritons created at the pumping spot (a) 
propagate along the channel (b) and accumulate on its right 
hand side (c). The fc-space behavior shows the decrease of 
the conensate energy (b,d,f). 



Gross-Pitaevskii treatment was developed where Boltz- 
mann scattering rates are dynamically calculated from 
the mean- field wavefunctions. The theory was successful 
in the description of micropillars with a small number of 
confined states, but it is not obvious how efficient the the- 
ory would be in the description of many different modes 
in a large reciprocal space. 

Comparing our results to the Monte Carlo approach 
developed in Ref. [sl], it should be noted that we are 
seeking to address a different problem. In Ref. [s^ the 
focus was on how quantum fluctuations affect the co- 
herence properties of parametric oscillators. Here the 
stochastic clement describes fluctuations due to phonons 
and we are interested in how they can cause relaxation of 
propagating condensates. To address coherence proper- 
ties, an approach based on a matrix of correlators would 
be more appropriate [s^ li^l. Unfortunately, such an 
approach is very demanding computationally, requiring 
matrices of size N x N, where N represents the number 
of states in reciprocal space. 

Finally, while we have considered the energy relaxation 
of exciton-polaritons via phonon interactions, wc should 
note that another relaxation mechanism is provided by 
scatterin g pr ocesses involving hot excitons with large mo- 
mentum |18j . Such hot excitons are typicall y cr eated in 
non- resonant /incoherently pumped systems |36| but can 
be neglected under resonant coherent excitation which 
we consider in this manuscript. In principle, the descrip- 
tion of scattering with high momentum excitons can be 
accomodated within our formalism. This would allow the 
theoretical study of the interplay between both exciton 
mediated and phonon mediated scattering processes in 
extended systems and is an important direction for fu- 
ture research. 

Conclusion. — We have derived a stochastic Gross- 
Pitaevskii equation, where the energy relaxation of 
bosons is provided by coupling to an incoherent field, 
treated as stochastic variable. As an example, we ap- 
plied the theory to the modeling of exciton-polaritons 
in a ID microwire with a potential gradient. The partial 
thermalization of exciton-polaritons is observed, together 
with their trapping in the real space. This result is of 
particular relevance to a variety of recent experiments in 
exciton and exciton-polariton systems where the energy 
relaxation of propagating Bose-Einstein condensates was 
reported. 

The work was supported by Rannis "Center of Ex- 
cellence in Polaritonics" and FP7 IRSES POLAPHEN 
grants. I.G. Savenko acknowledges support of the Eim- 
skip foundation. 
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SUPPLEMENTARY MATERIAL 

In the present supplemental appendix we derive the stochastic Gross-Pitaevskii-type equation, which accounts 
for the phonon-mcdiated relaxation. In the paper we focus on polariton propagation in a ID wire provided by 
confinement in the lateral direction, thus the appendix also corresponds to the ID case. It should be noted, however, 
that generalization of the formalism for the 2D case and for any system of interacting bosons is straightforward. We 
introduce the quantum field operators for polaritons, connected with annihilation operators in reciprocal {k-) 
space by the Fourier transforms (J^ in short notation). 



n^.] = 4=T.^:^e-^'''; (8) 



' y/N ^ 

i'x = ^-i[afc] = ^^a,e^'=% (9) 

^ k 

where N is the discretization length. 

The Hamiltonian of the system can be represented as a sum of two parts, 

H^ni + n2- (10) 

The first term Hi corresponds to coherent processes, 

Hi = ^'^"I^fe + E {y-'^i^- + c^-^i^i^x^x) + p.e^^'^-^-"-*) (i'l + , (11) 

k X X 

where is the (non-parabolic) polariton dispersion, T4 is the potential profile in x-space, a is a constant describing 
the strength of polariton-polariton interactions. The last sum in the equation is the coherent pumping term. Here 
Vx is the intensity of pump, kp is the wavevector of the pump, corresponding to the inclination of an incident laser 
beam, and fvjjp is the pumping energy. 

To model the interaction with acoustic phonons, we introduce the Frohlich-type Hamiltonian l^, 31[: 



q q,k 

where the phonons described by operators fet, h^, unlike polaritons, are considered to be three dimensional. The 
phonon wavevector \s q = BxQx + <^y<lx + Gz<1zi where Cx, &y and are unit vectors: e.x is in the wire direction, 
is in the structure growth direction, and e.y is perpendicular to ex and ez- The phonon dispersion relation, 

hu^ = huJq^ + 9^ + 92, is determined by the speed of sound, u. is the exciton-phonon interaction strength. 



which can be calculated using the standard electron - deformation potential interaction |33l - l35l | . As a result of elastic 
scattering of a crystall, there appears an area with relatively higher density and thus higher polarizability. This area 
attracts an electron efficiently acting as a potential well. The relative change in volume V can me expressed through 
the vector of displacement Udir) of a point called the deformation vector: 

SV 

— =V.^,. (13) 

Then, the expression for the electron energy can be written in form: 

E'[ud{r^] = - 4V • Ud{r1, (14) 

where Eq is the energy of electron in a not deformed crystall; the parameter is the constant of the deformation 
potential for an electron. For the lattice with one atom in an elementary cell, the vector of deformation can be written 
in form 1351: 
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where C is a real unit vector; p is the mean density of the semiconductor material. Putting p4|) in (|15|) . we obtain 
the scattering strength of an electron on a phonon: 



Gq — ~*C?e 



^ m 

2pVu 



2pVu^ 



(16) 



The same but opposite in sign expression is true for a hole. Combining them and introducing the overlap integrals, 
we finally come to the expression: 



del\iiqx,qv)lliqz) - dhlll{qx,qy)lliqz) 



(17) 



where dh are the deformation potentials of the lattice induced by phonons at the points of location holes. The 
integrals I^^^'^\qx: , qy) and l'^^\qz) are the overlap integrals of the phonon wavefunctions with the electron and hole 
wavefunctions, respectively, in the in- plane and growth directions and can be evaluated following Ref. (3l|. [BSj: 

-3/2 

e{h), ■ ' 



^11 iqx,qy) 



lf\qz) 



1 



nie + nih 
7r2 



ql + q^as 



2 



■sm 



qzLz 



(18) 
(19) 



It should be noted that a peculiarity of exciton-polariton systems is the variation of the excitonic fraction of polaritons 
with the in-plane wavevector, k, given by the /c-dependent Hopfield coefficient. Such a dependence would introduce 
a fc-dependence of the non-linear polariton-polariton interaction term and polariton-phonon interaction term. So as 
not to make our theory too specific to the exciton-polariton case, we neglect such a dependence. This is nevertheless 
valid in exciton-polariton systems in cases where the excited polaritons have similar in-plane wavevectors. In the 
excitation scheme that we consider, exciton-polaritons always have a small in-plane wavevector, due to the energy 
relaxation that limits their kinetic energy. Under such an assumption, the difference in a theory of excitons and 
of exciton-polaritons lies in the calculation of the non-linear interaction strength, a, and of the phonon scattering 
coefficient, G^. The nonlinear constant a should be multipled by and G by X, where X 1/2 for low energy 
exciton-polaritons and 1 for excitons. 

In cases where one wishes to describe exciton-polaritons with very different in-plane wavevectors, a theory could 
be based on writing independent equations for the evolution of coupled excitons and photons. In such a case, the 
interaction of excitons with phonons is of the same form as written in the current theory. 

The polariton field dynamics is given by: 



dt 



where the effect of the acoustic phonons is given by the term 



(20) 



(21) 



The evolution of the acoustic phonons is determined by the equation 



dh^ i 
~di ~ h 



^2, 



(22) 



The formal solution of Eq. 



reads 



h^{t) = b^{Q)e-'^st _ * / c-'^ii'-'')AqAt')dt', 
" Jo 



(23) 
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where 



A.W = E4'+g. W it). 

k' 

Substitution of this expression into Eq. (|21l) gives: 



(G,4_,,,,^,,jO)e— +G;4(0)e-^ 



(24) 



(25) 



^^,^|G5H'e-*«-"|^ e-*'*^f(*-*'U5,(t')di' - ^ e*"f(*-*'U,Ji')rfi'| , 



q 

where we noted that G-q^^q^^q^ = Gq^^q^^q. and uj^q^^q^^q^ = ^q^-qy-q^- Remembering that phonons represent an 
incoherent thermal reservoir, we can replace the terms b-qj.,qy,q^ (O)e^'"'''* and 6t(0)e'"'''* by stochastic classical variables 
b-q^^q^^q^{t) and b*^ respectively. This approximation is the analogue of the Markov approximation within 

the Langevin approach, when phonons are assumed to have a randomly varying phase. The stochastic variables are 
complex numbers with real and imaginary components drawn from a normal (Gaussian) distribution, normalized as 
follows, 



{bt{t)bq^{t')) ^ nq^6q-q^S{t ^ t'); 

{bS)h^{t')) ^ {b}{i)bt,{t')) 



(26) 
(27) 



where is the number of phonons in the state with wavevector q determined by the temperature of the system. The 
summation of the stochastic terms over qy and qz in Eq. ([25]) can be made by noting that the sum of two normally 
distributed stochastic variables again gives a normally distributed stochastic variable. 



G-q:,.qy,q,b-q^^q^^qXt) — 71g^(i); 



x.qy.q^ q^o.qy 



,qAt) = r:M 



(28) 
(29) 



qy,qz 



where Tq^ and T*^ represent the temperature dependent or thermal part of polariton scattering on phonons. They 
are defined by the correlators: 



qm,qy,qi^qx,q'^^i^ ^ )i 



a:(i)r,i(0)= E \Gq.,qy,q. 

qy.qy 

{TqAt)Tq>At'))^(Tq:it)Tqlit'))=0. 

The time integrals in Eq. (|25p can be simplified using the fact that the function 



(30) 
(31) 

(32) 



^ ^'feif / / \G{q)\hm[uj{q)t]dqydqz 



entering Eq. ([20]) . is approximately independent of q^ in the range of Qx G (—10®, 10®) m ^, and thus ICq^ (t) w ICo{t). 

Further, within the Monte Carlo approach, the field operator can be treated as a classical variable for condensed 
polaritons, ipx = (^'x) (with the Fourier image ipk)- Then, physical observables are calculated over multiple realizations 
of the evolution dynamics with stochastic variable bk{t). We make a last notation 



Skit) - ^^+9=" (*) (/ A. (t'w - t')dt'^ , 



and finally obtain the stochastic Gross-Pitaevskii - type equation 



ih- 



dt 



: [EkiJk + Skit)] + 

■Y.{T-kit) + Tkrit)}c 



Vx - + a {ipxl^ 



i'x + p^e^'^^'^e-'""* 



(33) 



(34) 



'ikx 
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where we introduced the decay term —ihrftjjx/2 to account for the radiative decay of particles [19|. The function Sk{t) 
is identified as the stimulated part of the phonon-mediated relaxation, where the convolution integral is responsible 
for energy conservation during the density-stimulated phonon-assisted process of energy relaxation. 

Equation p4p is the main result of the Letter. It represents the dissipative GP equation and can be applied to 
various bosonic systems. As an example of application, here we present the figure which illustrates energy relaxation 
in a parabolic trap (see Fig. [3]). 




-60 -40 -20 20 40 60 





FIG. 3: Relaxation of energy of exciton-polaritons in the parabolic trap in the regime of cw excitation in a quantum wire due 
to phonon-assisted processes in x- and fc-space for the time range: 0-50 ps. White curves show the Harmonic-oscillator like 
potential profile in the real space and free exciton-polariton dispersion in the reciprocal space used in calculations. Coherent 
pumping represents a Gaussian in x-space centered around x = /im. The inclination angle of the pump in the units of 
wavevector is fcp « —0.9 fim~^, the energy hwp = 0.5 meV. Polaritons created at the pumping spot, first, propagate until they 
meet the potential profile and then along the potential gradient and accumulate in its bottom (a). The fc-space behavior shows 
the decrease of the conensate energy (b). 
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